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Abstract. The complex band structures calculated using the Extended Plane Wave 

Expansion (EPWE) reveal the presence of evanescent modes in periodic systems, never 
Q ' predicted by the classical uj{k) methods, providing novel interpretations of several 

O ^ ' phenomena as well as a complete picture of the system. In this work we theoretically 

and experimentally observe that in the ranges of frequencies where a deaf band is 
\l ' traditionally predicted, an evanescent mode with the excitable symmetry appears 

changing drastically the interpretation of the transmission properties. On the other 
^o I hand, the simplicity of the sonic crystals in which only the longitudinal polarization can 

VO ■ be excited, is used to interpret, without loss of generality, the level repulsion between 

symmetric and antisymmetric bands in sonic crystals as the presence of an evanescent 



C^^ ' mode connecting both repelled bands. These evanescent modes, obtained using EPWE, 



explain both the attenuation produced in this range of frequencies and the transfer of 
symmetry from one band to the other in good agreement with both experimental results 
and multiple scattering predictions. Thus, the evanescent properties of the periodic 
system have been revealed necessary for the design of new acoustic and electromagnetic 
applications based on periodicity. 
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1. Introduction 

Since tlie pioneering works of Yablonovitch [H |2] and John [3J simultaneously discovered 
the possibihties to control the flow of light in periodic distribution of dielectric materials, 
a drastic increasing interest appeared in the analogous structures to control both the 
elastic and acoustic waves using the well-known Phononic [U [5] or Sonic Crystals 
O [HI m [To] respectively. Phononic Crystals are periodic distributions of two elastic 
materials, however in the case that one of these materials is a fluid the system is called 
Sonic Crystal (SC). In the past twenty years the exploitation of the particular dispersion 
relation of these structures has revealed very interesting physical properties [H [5], [6l [TJ [HJ 
[a [101 [III [121 [131 [11 [ig [ig [ni [111 [191 [201 [211 [221 [231 ^ 

The existence of ranges of frequencies, called Bandgaps (BGs), in which non propagating 
modes can be excited in the system, has been observed in a wide range of frequencies 
due to the scalability of the systems [131 HH [2H]- The presence of the EG has motivated 
the analysis of both the localized states in point defects [HI [121 113 [13 [22] and the 
guidance of waves thorough a linear defect constituting a waveguide in the periodic 
system [TH [HI [25]. However, other properties of the dispersion relation have been used 
to control the wave propagation through these periodic structures. On one hand, the 
linearity of the dispersion curves at subwavelength regimes has motivated the application 
of homogenization theories in order to deflne effective parameters of the system [161 El] • 
On the other hand, the curvature of the isofrequency contours has been used to control 
the spatial dispersion of the waves inside the periodic structures in order to obtain both 
self-collimating effect [261 [2Z] and negative refraction [2Ql [21] . 

Recently, the possibility to control the evanescent properties in periodic composites 
has shown several interesting possibilities in both photonic [351 [36] and phononic 
[371 [381 [39] crystals. Imaging with super resolution, [lOl [H] that is, overcoming the 
diffraction limit can be obtained by restoring all the evanescent components of a near- 
fleld image. This can be achieved by the coupling of the evanescent modes with 
other mechanism leading to their ampliflcation in order to successfully transport the 
information carried by the evanescent waves through the system. In this way, the 
evanescent properties of the periodic system have been revealed necessary for the design 
of new acoustic and electromagnetic applications. 

The dispersion relation of these periodic structures has been traditionally obtained 
using u{k) methods, for instance, the plane wave expansion (PWE) [7]. Using the 
PWE method, the BGs can be deflned as ranges of frequencies in which no modes are 
excited inside the periodic structure. However, BGs should be deflned as ranges of 
frequencies where only evanescent modes can be excited [121 HSj- Recent works have 
shown in both regimes, acoustic [31] and elastic [H], that an extension of the PWE 
method, called the Extended Plane Wave Expansion (EPWE), is an efficient tool for 
analyzing the evanescent properties of the periodic systems. This methodology inverts 
the u{k) problem to the k{u) problem by fixing both the angular frequency u and 
the direction of incidence a {k = ka), k being the eigenvalue of the problem. The 
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supercell approximation has been introduced in the EPWE to theoretically characterize 
the evanescence of both single |30] and double |32] point defects. The evanescent decay 
predicted by EPWE for both complete and defective SC is in very good agreement with 
the experimental measurements [SQl |3T1 [32]. The picture of the propagation properties 
of the periodic systems provided by the EPWE method improves and completes that 
given by the PWE. 

The goals of the paper are two-folds. First we observe the evidence of 
evanescent modes in a sonic crystal waveguide (SCW) using the EPWE with supercell 
approximation. Second, we analyze the evanescent properties of stubbed SCW. The 
presence of the stub, breaks the symmetry of the system and fundamental phenomena 
can be studied. The antisymmetry of the system introduces several splittings in the 
band structure characterized by an interchange of the properties between the repelled 
bands. In this work we study the nature of the level repulsion in two-dimensional 
periodic systems making use of evanescent waves. 

We find the evanescent modes of fundamental interest for the correct understanding 
of the wave propagation properties in both linear waveguides and in general structures 
based on periodicity. We use a SCW because it allows us to obtain high precision 
measurements inside the periodic structure where the evanescent properties are excited. 
As it has been previously shown [31], the interest of SCW stands on the effects due to 
the periodicity that can not exist in rigid walls waveguides. Due to the linear defect 
in the SCW interesting physics appear in the BG. In this work we theoretically and 
experimentally observe that for the deaf bands [HI Il5l HE], where the real part of 
the complex band structure represents antisymmetric modes, evanescent modes with 
the excitable symmetry can appear in the waveguide. In contrast to the classical 
interpretation of the attenuation in a deaf band, based on geometric arguments, in 
this work we find that the attenuation is governed by a decay rate related with the 
imaginary part of the complex band structures obtained using EPWE. Independently, 
we have obtained the same value of the decay rate using other methodologies, as the 
Multiple Scattering Theory (MST) ^7\ SHI HHl EQ], in very good agreement with the 
experiments. 

On the other hand, the mode conversion in phononic crystals thin slab has been 
recently observed [51]. This conversion depends on the symmetry of the system with 
respect to the direction of the incident wave. In these systems, the excited modes 
either couple or not with other ones, being the coupling observed by both a splitting 
between bands in the band structure and a transfer of the symmetry from one band to 
the other. In the case in which no coupling occurs, the bands simply crosses between 
them. Bavencoffe [52] established a link between the attenuation of the ultrasonic wave 
observed in the case of a limited grating and the values of the imaginary part of the 
wave number in a stop band computed for an infinite grating. In other works, Achaoui 
et al. [53] , have recently observed in phononic crystals with some freedom of anisotropy 
that when a band mostly polarized in-plane is close to a band mostly polarized out-of- 
plane, a phenomenon of repelling can occur and in some instances it introduces a local 
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band gap. Moreover, this interaction is accompanied by a transfer of the polarization 
state from one band to the other. These levels repulsion which avoids crossing in the 
distribution of eigenvalues is well known within the physics community [Ml [551 l56] . 
However, there are only few works in the literature that analyzes and discusses this 
phenomenon in two-dimensional phononic crystals [571 EH l53] . 

The simplicity of the SC in which only the longitudinal polarization can be excited, 
is used to observe without loss of generality the level repulsion between symmetric 
and antisymmetric bands. In this work we show that the level repulsion between the 
symmetric and antisymmetric modes is the result of the presence of an evanescent 
mode. These evanescent modes explain both the attenuation produced in this range 
of frequencies and the transfer of symmetry from one band to the other. 

The work is organized as follows. First of all both the theoretical techniques and the 
experimental set up are shown in Section O In Section |3] we present the classification of 
the several evanescent modes predicted by the EPWE showing the main characteristics 
of the complex band structures of SCW. The evanescent behaviour of modes in SCW 
is also analyzed both analytically using MST and experimentally. Section Hj analyses 
the stubbed SCW. The evanescent bands, non predicted by the traditional PWE, are 
very relevant in these systems. The interaction of the evanescent modes with cavity 
modes in the waveguide is shown in Section 14.11 In Section 14.21 the stubbed SCW are 
used to show the effect of the level repulsion in periodic systems, showing both the local 
bandgaps and the symmetry transfer between symmetric and antisymmetric bands due 
to the presence of evanescent modes. The concluding remarks of the work will be shown 
in Section [5l 

2. Theoretical techniques and experimental set up 

The propagation properties of periodic materials can be analyzed solving for both the 
eigenvalue and the scattering problems. The former can be used to obtain the dispersion 
relation of a system whereas the latter can be used to obtain the scattering of waves in 
finite structures. Both methodologies are independently used in this work to analyze the 
properties of an acoustic waveguide. The eigenvalue problem k{u) is solved in this work 
using the EPWE with supercell approximation [32]. To compare the results with the 
classical uj{k) some calculations using the PWE method [7] will be shown. To solve the 
scattering problem we use the MST [171 HH HHl [50], which is a self-consistent method 
used to analyze the scattered field of an arrangement of scatterers considering all the 
orders of scattering. 

The experiments showed in this work have been performed in an echo-free chamber 
sized 8x6x3 m^. To obtain the experimental dependence of the pressure all along the 
SCW, we measured the pressure field inside the guide. We build a finite SC bunging 
the rigid scatterers on a periodic frame (Figure [I]) and the SCW is generated removing 
the central row of the complete structure. We notice that this system avoids the excess 
attenuation effect [33]. The finite 2D SC is made of 7x7 rigid cylinders arranged as a 
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square array of lattice constant a. The radius of the cyhnders used in the experiments 
is r = 9 cm, and the lattice constant of the SC is a = 19 cm. With these parameters, 

2 

the filling fraction of the finite SC is // = ^ ~ 0.71. The dimensions of the system are 
large enough for the microphone to be placed inside the waveguide. The microphone 
used is a prepolarized free-field 1/2" Type 4189 B&K. The diameter of the microphone 
is 1.32 cm, which is approximately 0.07a, therefore a low level of influence over the 
pressure field measured is expected. The source is placed 1 m away from the SCW 
launching white noise. 

In this work a 3D computer-controlled automatic positioning system together 
with an automatized acquisition system, called 3DReAMS (3D Robotized e-Acoustic 
Measurement System), has been used to obtain the pressure field inside the waveguide. 
The 3DReAMS system is capable of sweeping the microphone through a 3D grid of 
measuring points located at any trajectory inside the echo-free chamber. The robot 
was controlled by the motion controller card of National Instruments PCI 7334. We 
analyzed the acoustic field inside the guide of both complete SCW and stubbed SCW 
by moving the microphone in steps of 1 cm. 




Figure 1. (Color online) Experimental set up. Picture of the frame in which the SCW 
has been hung and the complete structure with 7x7 rows in which the central one is 
removed in order to produce the waveguide. 



3. Propagating and evanescent modes in Sonic Crystals Waveguides 



Consider a periodic arrangement of rigid scatterers embedded in air {pair = 1-29 kg/m^, 
c=340 m/s) with square periodicity which lattice constant is a and the filling fraction 
// = nr'^/a'^ = 0.717. All the physical parameters presented in the work will be shown in 
reduced magnitudes: reduced frequencies (\1/ = ua/cair, v is the frequency) and reduced 
Bloch vector [K = ka/ {27r)). The band structure for this periodic array is shown in 
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Figure |2^. The bandgap of the complete SC (without hnear defect) is ranged between 
\E'i = 0.33 and \l/2 = 0.8. In this Section we analyze the evanescent and propagating 
properties of a SCW generated by removing the central row of the complete structure. 
EPWE shows novel properties of the guided and non-guided modes in the waveguide 
that would not be observed otherwise using PWE. 

3.1. Preliminary analysis: Transmission properties of a Sonic Crystal Waveguide. 
Motivation 

The generated SCW is completely equivalent to that used in a previous work [58] . 
Figure [2]d shows the band structure along the FX direction for the SCW for these 
frequencies. The band structure has been calculated using the PWE with the supercell 
approximation. 1105 plane waves have been used for the calculations. The supercell is 
shown in the inset of Figure |2}3. 

Because of the linear defect, the band structure in Figure Eb shows five guiding 
bands in the bandgap. The red continuous and blue dashed lines in Figure [2)3 correspond 
to localized modes in the straight waveguide. The Fourier transform of the eigenvectors 
yields the pressure field inside the SCW. Analyzing this field |58j, one can observe 
that the modes corresponding to the red continuous lines produce symmetric acoustic 
fields with respect to main direction of symmetry FX, whereas the acoustic fields 
corresponding to the blue dashed lines are antisymmetric. Thus the symmetric modes 
can be excited by a longitudinal incident wave and the antisymmetric modes will not 
contribute to the transmission because they cannot be excited. The later modes are 
usually called deaf bands and they have been analyzed in the literature P HHl HE] . 

From the previous discussion it is expected that between symmetric guiding bands 
should appear a dip of transmission (grey areas in Figure [2]d). In order to observe this 
property, one should additionally solve for the scattering problem. We analyze several 
SCW with different number of rows using MST. We evaluate the Insertion Loss (IL) 
[9], which corresponds to the attenuation spectrum evaluating the difference between 
the sound level recorded in the same point without and with the SCW. The comparison 
between the band structure and the IL spectrum shown in Figure [2b indicates that 
the symmetric bands (red continuous lines in Figure [2)3) contribute to the transmission 
whereas, as it was expected, between this symmetric bands, attenuation peaks appear 
in the IL spectrum because of the antisymmetric modes (blue dashed lines) cannot be 
excited. 

Figure [2t clearly shows the two attenuation peaks produced by the non excited 
antisymmetric bands in the waveguide for a high number of rows. However, in the case 
of low number of rows the attenuation level due to the first antisymmetric band is very 
low in comparison with the attenuation level of the second peak (red continuous line 
in Figure [2}:;) . Moreover, the attenuation level of the first attenuation peak increases 
with the number of rows whereas the second peak seems more stable. Thus several 
fundamental questions reach from this plot. First of all, if band structures predict stop 
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bands for these frequencies (grey areas), what kind of modes can propagate depending on 
the length of the waveguide? Why does the second antisymmetric mode present larger 
and more stable attenuation level than the first one? The answers of these questions are 
related to evanescent waves and cannot be obtained from the results of either the PWE 
nor the MST. A more accurate approach should be used to understand the physics of 
these systems. 
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Figure 2. (Color online) Study of the transmission properties of the SC and 
SCW. (a) Bands structure for the complete SC made of rigid scatterers embedded 
in air with square periodicity which lattice constant is a and the filling fraction is 
// = 7rr^/a^ = 0.717. Grey area represents the bandgap. (b) Bands structure of the 
SCW at rX direction calculated using PWE (1105 plane waves) for the frequencies 
in the bandgap of the complete structure. Grey areas represent non guiding bands. 
Red continuous (Blue dashed) lines represent symmetric (antisymmetric) bands, (c) 
IL prediction calculated using MST. The different lines represent the IL predictions 
for SCW with several number of rows (see legend). 



3.2. Complex band structures in Sonic Crystals Waveguides 



Brillouin [l2] revealed that waves in the bandgap present evanescent behaviour, i.e., 
they should be characterized by a complex value of the wavevector in such a way the 
amplitude of the mode should be modulated by an exponential with negative exponent. 
Therefore the mode should decay as it penetrates in the crystal. The elegant and 
intuitive explanation of Joannopoulos |13] revealed that the decay rate should be greater 
for frequencies closer to the center of the BG. Classical methods do not show this kind 
of properties because they interpret the BG as ranges of frequencies where no modes 
can be excited. 

The Joannopoulos predictions have been recently observed both theoretically 
and experimentally [3T| [30] for SCs. Motivated by the properties of the complex 



Evanescent modes and level repulsion in SCW 



(b) 
l!::::;:ttww» 

J,J.„......— """ 

9!:::II,|l;;j::...... 




0,2 

lm(K) 



Figure 3. (Color online) Analysis of the band structure of the SCW. Grey area 
represents the bandgap. Reduced magnitudes have been used in the plots, (a) 
Calculation using EPWE. Red triangles represent modes with Im(K) ^ and 
Re{K) = 0, black circles represent modes with Re{K) ^ and Im{K) = 0, green 
squares show modes with Re{K) — 0.5 and Im{K) ^ and blue diamonds represent 
modes with < Re{K) < 0.5 and Im{K) ^ 0. (b) Bands structure calculated 
using PWE. Inside the bandgap, open circles (black circles) represent antisymmetric 
(symmetric) modes. The supercell used in both cases are shown in the insets of (a) 
and (b). 



band structures and by the interesting questions related to evanescent modes presented 
in the previous Section, we study the evanescent and guiding properties of SCW. The 
conclusions of this work are general and can be applied to other kind of periodic systems. 
Due to the freedom of the eigenvalues in the EPWE method (complex valued k), 
a more general band structure than the classical ones are obtained [311 El- I^i this 
Section we classify the propagating and evanescent modes appearing in the complex 
band structure resulting from the EPWE method. The classification shown in this work, 
inspired by the works of Bavencoffe et al. [521 [59] , helps to a better understanding of 
the complex band structures. The modes are classified following the next restrictions: 
(?) The classical band structures correspond to modes characterized by values of Re{K) 
in the Brillouin zone and Im{K) = 0. Modes with these properties are shown in 
this work with black filled circles, (u) The modes characterized by Im{K) > and 
Re{K) = are shown with red filled triangles. These modes represent connections 
between propagating bands in the F point. (Hi) The modes characterized by Im{K) > 
and Re{K) = 1/2 (respectively, Re{K) = 1/a/2) are shown with green filled squares. 
These modes represent connections between propagating bands in the X (respectively, 
M) point, (iv) The modes with Re{K) in the first Brillouin zone but with Im{K) ^ 
are shown in blue filled diamonds. These modes belong to evanescent connecting bands 
between bands with the same symmetry of the acoustic field. For this last group we 
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constrain the plots for modes with Im{K) < 0.4. 

Following the previous classification we have plotted the complex band structure 
calculated using the EPWE with supercell approximation for the SCW we are dealing 
with. The convergence of all the calculations have been carefully analyzed selecting the 
adequate number of plane waves for each case. 1369 plane waves have been used for the 
calculation. Figure [3^ shows the complex band structure for the SCW at FX direction. 
The values of the real part of the wave vector (central panel of Figure [3^) are restricted 
to the first Brillouin zone, whereas the imaginary part (right and left panels of Figures 
[3^) is arbitrarily displayed. It is worth noting that due to the frequency discretization 
to solve the EPWE, it is difficult to find fiat bands in the band structure. Thus, the 
frequency step used in this work has been specifically selected for each calculation. The 
grey area shows the range of frequencies of the bandgap of the perfect complete SC. 
The supercell used for the calculations is shown in the inset of Figure [3ti. 
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Figure 4. (Color online) Analysis of the propagating and evanescent modes of the 
SCW predicted by EPWE (modes at * = 0.484 *2 = 0.619 and ^3 = 0.731). (a) 
Complex bands in the band gap region, (bl)-(dl) Modulus of the acoustic field for the 
real eigenvalues. (b2)-(d2) Modulus of the acoustic field for the complex eigenvalues. 



Comparing the band structures calculated using EPWE and PWE in Figures [3^ and 
|3b respectively, one can observe that EPWE reproduces the same guided modes in the 
bandgap as in the case of PWE. However two important characteristics should be noted. 
First, the complex band structure shows that two consecutive symmetric (respectively, 
antisymmetric) real bands are connected by a symmetric (respectively, antisymmetric) 
imaginary band indicated with black continuous (respectively, dotted) arrow in Figure 
|3K. Therefore the complex connections preserve the symmetry of the bands which are 
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connecting. Second, as a consequence of the connections and due to the folding of the 
bands, an additional band not predicted by PWE appears in the bandgap (see real and 
imaginary bands represented with blue diamonds). Let us analyze in detail these two 
features in the next lines. 

At r or X points where the folding of the bands is expected, an evanescent 
connection between the bands should appear to preserve the overall number of bands 
at any frequency. In order to observe this, we fix our attention in the bandgap for the 
complete periodic system. Figure HIi shows the complex band structure for this range 
of frequencies. Let us analyze the first symmetric (respectively, antisymmetric) guided 
band labelled as srl (arl) in Figure Hti. Note that this band is connected to the second 
symmetric (antisymmetric) band, sr2 (respectively, ar2), by means of a complex band 
with non zero imaginary part sil (ail). It is worth noting that both the maximum value 
of the band srl {arl ) and the minimum value of the sr2 {ar2) are at point X. Thus the 
complex band that connects these two symmetric (antisymmetric) real bands presents 
a constant value of the real part whereas the imaginary part follows a connecting path 
between the two real bands preserving the symmetry of the band. 

In the range of frequencies between the two symmetric (antisymmetric) bands, PWE 
predicts that due to the antisymmetry (symmetry) of the acoustic field, stop (guiding) 
bands are expected. Here, between the symmetric bands srl and sr2, an evanescent 
symmetric bands, sil, is excited whereas between the two antisymmetric bands arl and 
ar2, a symmetric real band is excited, sr2. Thus, the interpretation of EPWE gives 
a complete picture of the physics of the problem. The stop bands can be interpreted 
as ranges of frequencies where symmetric evanescent modes are excited; the decay rate 
of the mode inside the SCW being linked to the imaginary part of the complex wave 
vector. 

To interpret these results we have selected two different modes at two different 
frequencies, one at ^i = 0.484 {Ki = 1/2 + 0.049«) and the other one at ^2 = 0.619 
{K2 = 0.35). The analyzed modes are marked in the band structure in Figure H^ 
with red squares. For the analysis we only select the first complex band because 
the complex bands with high values of Im{K) do not contribute to the transmission 
properties of the SCW. Figures Hbl and|4b2 (jlbl and|4l:;2) show the acoustic field of the 
mode for the real and complex eigenvalue at the frequency \l'i {"^2) respectively. The 
acoustic field has been obtained from the Fourier transform of the eigenvectors for the 
corresponding frequency. We would like to note that the real mode at the band arl (sr2) 
is exactly the same as the predicted using PWE with a value oi K = 0.33 {K2 = 0.35) 
presenting antisymmetric (symmetric) acoustic field. Oppositely, the complex mode 
only predicted by EPWE at the band sil (ail) is characterized by a value of the 
Bloch vector Ki = 1/2 + ^0.049 {K = 1/2 + ^0.15) with a symmetric (antisymmetric) 
evanescent acoustic field. Thus, because of the symmetry of the SCW, only the modes 
with symmetric acoustic field are excited. Therefore for the frequency \l/i the stop band 
is characterized by the evanescent mode shown in Figure Ilb2 whereas at frequency \l/2 
the symmetric propagating mode shown in Figure Htl is expected. 
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Figure 5. (Color online) Multiple scattering analysis of the excited modes in a 
waveguide made by removing one row of scatterers located at y/a = in a SC 
constituted of 10x10 cylinders. (al)-(a2), (bl)-(b2), (cl)-(c2) show the acoustic field 
inside the waveguide and the longitudinal cut at y = (black continuous line) for 
the frequencies 0.484, 0.619 and 0.713 respectively. Red line shows the exponential 
fit of the red circles that represents the decay of the mode. Green line represents 
the decay predicted using the Im{K) from EPWE. Open black circles represent the 
measurements in the experimental set up (see Section [2]). 



The connection between the bands sr2 and sr3 is different than in the previous 
cases. The maximum value of the band sr2 is at F point whereas the minimum value 
of the band sr3 is at X point. Then, the evanescent connection between these bands 
presents both a real part that crosses the path FX, sre, and an imaginary part, sie, 
connecting the bands in the imaginary part. This connection produces an additional 
band that has never been predicted by PWE. To analyze this situation we have selected 
two modes at frequency \E'3 = 0.731, one with real value of the Bloch vector, K = 0.44, 
at band ar2 and another one with complex value of the Bloch vector K^ = 0.33 + «0.21, 
being the real part of K in the band sre and the imaginary part in the band sie. The 
acoustic fields of the real and the complex modes are shown in Figure HJil and |lli2 
respectively. One can observe that the acoustic field corresponding to the real mode, 
shown in Figure HJil, is antisymmetric and cannot be excited in the SCW. However the 
evanescent field produced by the complex mode. Figure 11112, shows a symmetric acoustic 
field and, as a consequence, it is an excitable mode in the SCW. 

The EPWE predictions allow us to affirm that the two stop bands in the SCW (grey 
areas in Figure [2)d) are characterized by two evanescent symmetric modes. The first stop 
band is characterized by an evanescent band that presents lower values of the imaginary 
part than in the case of the second stop band. This feature explains both the dependence 
of the IL spectra on the number of rows and the stability of the attenuation level of 
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the second peak shown in Figure [2]3. However, due to the assumption of periodicity, 
the EPWE predictions are for an infinite SCW. In order to both check and clarify the 
previous results we have calculated the acoustic fields in a real waveguide made from 
a complete SC of 11 x 11 rows in which we have created a linear defect in the central 
column (columns are defined in the direction of incidence). 

Figure [5] shows the MST calculation for the three frequencies previously analyzed: 
(a) ^1 = 0.484, (b) ^2 = 0.619 and (c) ^3 = 0.731. Figures marked with I's shows 
the absolute value of the acoustic pressure field inside the waveguide, whereas the 
Figures labelled with 2's show the longitudinal cut of the field inside the waveguide 
for y/a = 0. The MST calculation predicts acoustic fields inside the waveguide similar 
to the symmetric ones obtained from the Fourier transform of the eigenvectors of the 
EPWE for the analyzed frequencies. This means that the excited acoustic field in the 
SCW corresponds to the first symmetric mode predicted by EPWE. Thus, the mode at 
\1' =0.619 should propagate through the waveguide whereas the modes at \l/ =0.484 and 
\1' =0.731 should be evanescent and the decay should be comparable with the value of 
the imaginary part of the K. 

To analyze both the propagating and the evanescent behaviour of modes we study 
the acoustic field in the longitudinal cut at y/a = 0. Figure [5]d2 represents the cut 
corresponding to the propagating mode at \E' = 0.619. One can observe that the field 
does not present any decay and the wave is guided through the waveguide. Oppositely, 
Figures |5^2 and |5t2 represent the cuts for the evanescent modes at 0.484 and 0.731. 
These cuts show clearly the decay produced because of the evanescent nature of the 
excited mode. The experiments (open black circles) are in good agreement with the 
evanescent behaviour of the modes. 

In order to check that effectively the excited modes in the stop bands corresponds to 
the evanescent mode obtained using EPWE, we have fitted the decay of the modes to an 
exponential ae~*^. The fit of the mode at \1' = 0.484 using the maximum points shown 
in Figure [5^2 with red circles, gives a value of the decay rate Im{K)fit = 0.053 ± 0.004. 
This value is in complete agreement with that predicted by EPWE, Iin{K) = 0.049. 
On the other hand, analogous procedure gives a decay rate for the mode at \I' = 0.731 
equal to Im{K)fit = 0.21 ± 0.03 in complete agreement with Im{K) = 0.21 (Figure 
|5fc2) predicted by EPWE. 

4. Evanescent coupling and level repulsion states in SCW 

The interaction between stubs or point defects with waveguides has been extensively 
analyzed in the literature |58l [34l |60]. The presence of the stub or a defect in a SCW 
alters drastically the transmission spectrum. Although the general characteristics of the 
transmission properties of the SCW are preserved once the stub or the defect is created 
in the waveguide, narrow attenuation peaks appear in the guided bands occurring at 
the resonance modes of the stub. The stub modes are evanescent and strongly localized 
by the surrounding periodic structure. 
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Figure 6. (Color online) Theoretical and experimental analysis of the stubbed 
waveguide using EPWE and MST. (a) IL spectrum evaluated at the end of the 
waveguide. In the central panel blue continuous line (red dashed line) represent 
the MST predictions and in the left panel blue dots (open red circles) represent the 
experimental results of the SCW (stubbed SCW). Peaks labelled as R (LR) represent 
the effect of the resonances of the stub (effect of the level repulsion between bands 
with different symmetry). The right panel represents the real part of the complex 
band structure for the SCW. (b)-(f) show the details of both the hybridization of the 
resonances and the guide or evanescent modes and the repulsion between bands (see 
text). 



In this Section, we investigate both theoretically and experimentally the occurrence 
of the effects produced by a stub in a SCW. Stubs are generated in this work by removing 
a scatterer in the side-wall of the waveguide. We observe two different effects due to 
the stub. On one hand, EPWE predicts attenuation dips due to the interaction of 
both the guided and the evanescent modes with the resonances of the stub. This result 
constitutes a proof of the importance of the k{u) methods because they predict the 
evanescent bands and, as a consequence, the interaction between the evanescent modes 
and the resonances of the stub. The attenuation peaks in Figure [6^ marked with R 
(resonance) correspond to this first effect and they will be discussed in detail later. 

On the other hand, the presence of the stub changes the symmetry of the system, 
therefore antisymmetric modes can be excited even by symmetric longitudinal incident 
waves. This produces the appearance of repulsion states between symmetric and 
antisymmetric bands. In this work we observe that due to the local breaking of 
the symmetry because of the presence of the stub, a repelling state appears between 
symmetric and antisymmetric bands accompanied by a transfer of the symmetry of the 
modes from one band to the other. Chen et al. showed that depending on the symmetry 
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of the system with respect to the incident direction of the incident wave, shear-horizontal 
modes either couple or not with the Lamb wave modes which are polarized in the sagittal 
plane. The coupling can be observed by a splitting between bands in the band structure 
and by a transfer of the symmetry from one band to the other. In the case in which no 
coupling occurs, the symmetric Lamb waves band simply crosses of the shear-horizontal 
band. Similar results have been recently observed in ID phononic crystals |52l|59], and 
anisotropic phononic crystals [53]. We observe in this work a new interpretation: The 
repelled bands are connected by an evanescent mode which is the responsible of the 
transfer of the symmetry from one band to another. The attenuation peaks in Figure 
Eh. marked with LR (level repulsion) correspond to this second effect and they will be 
discussed in detail later. 

The real part of complex band structures of the stubbed waveguide are also shown 
in the right panel of the Figure EK. One can observe the correspondence between the 
attenuation peaks, R and LR, predicted by MST (central panel of Figure [6^1) and the 
singular points of the band structure. Left panel of Figure Eh shows the experimental 
data (see Section [2] for the details of the experimental set up) in good agreement with 
the theoretical predictions. Let us study the nature of each of these attenuation peaks. 

4.I. Coupling between guided and evanescent modes with stub resonances 

The stub acts as a resonant cavity in the waveguide producing flat bands due to the 
localized mode in the stub. If this resonance occurs at the same frequency as a guiding 
or evanescent bands, then the hybridization of the two bands produces an attenuation 
peak in the spectrum. Resonances of the stub, calculated using EPWE at point X, can 
be observed in Figure [3 It is worth noting that the frequencies of these resonances 
coincide with two symmetric guiding bands and one evanescent band of the waveguide. 




Figure 7. (Color online) Absolute value of the acoustic field for the stub resonances, 
(a) Rl, 1'i==0.3812. (b) R2, *2=0.605. (c) R3, *3=0.72. 

Points Rl (^ = 0.3812) and R2 (^ = 0.605) in Figure Eh coincide with the 
hybridization of the guiding bands and the resonance bands of the stub (fiat band). A 
detailed image of the hybridization effect at these frequencies can be seen in Figures Eb 
and Eh- These hybridizations explain the attenuation peaks in the IL spectrum labelled 
as Rl and R2 in Figure Eh in good agreement with both the analytical predictions and 
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the experimental results. The shape of the attenuation peaks assembles the Fano-like 
resonances previously reported [61] • The slight disagreements between the experimental 
data and the theoretical predictions can be due to the sensitivity of the experimental 
data to the size of the stub. 



1t0.605 




Figure 8. (Color online) Coupling of guided and evanescent modes with stub 
resonances, (a) Coupling of the resonance at ^ = 0.605 with the guided mode at 
the same frequency, (b) Coupling of the resonance at ^ = 0.72 with the evanescent 
mode at the same frequency. 



The interaction between the guided mode in the SCW and the resonance of the 
stub at point R2 is analyzed in Figure |Ht. The calculations have been carried out using 
the MST. We study the propagation of waves at frequency "^ = 0.605 inside a SCW 
7a long with a stub. The position of the stub has been changed all along the SCW in 
order to observe that the guided mode excites the resonance of the stub independently 
of the distance between the site of the stub and the entrance of the SCW (left side of 
the SCW). The wave penetrates in the waveguide exciting the symmetric guided mode. 
Once the guided mode finds the stub, the resonance of the stub is excited, localizing the 
field in that region. One can compare the field around the stub in Figure [Hk, calculated 
using MST, with the acoustic field of the resonance of the stub in Figure [7)d calculated 
using EPWE. It is worth noting that the resonance of the stub is excited practically 
with the same intensity in all the locations of the stub due to the low losses of the guided 
mode. Analogous effect occurs for the point Rl (^ = 0.3812). 

Point R3 (\E' =0.72) corresponds to the interaction of an evanescent mode with the 
resonance of the stub at this frequency. The detailed band structure for this frequency 
is shown in Figure [6]f. Note that for this case no hybridization is produced and two 
bands can be excited. Once the wave penetrates in the SCW the evanescent mode is 
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excited and, if the stub is as close as the evanescent mode can travel through the SCW, 
the stub mode should be excited. Note that these results are not predicted by the 
classical methods, therefore they cannot be explained with the PWE, but EPWE gives 
the correct understanding of the phenomenon. 

MST simulations have been carried out in order to observe this coupling between 
the evanescent mode and the stub resonances at \1/ =0.72. Figure [8]d shows five stubbed 
waveguides with increasing distances between the sites of the stub and the entrance 
of the SCW. In all cases, the excited evanescent mode presents a symmetric acoustic 
field equivalent to that predicted by EPWE in Figure |lli2. Once this mode arrives to 
the stub, it excites the stub resonance. However it is worth noting that the mode is 
evanescent and, as a consequence, the intensity of the localized modes in the stub should 
decrease as the stub is far from the entrance of the SCW. The MST calculations predicts 
this behaviour. For the last case, the distance between the stub and the entrance of 
the SCW is larger than the penetration distance of the evanescent mode and, as a 
consequence, any stub resonance can be excited. 

4-2. Level repulsion states 

In addition to the attenuation peaks due to the resonances, the IL spectrum of the 
stubbed waveguide (Figure E^i) predicts two additional attenuation peaks referred to as 
LRl (\1/ =0.429) and LR2 (\1/ =0.544). These peaks have different phenomenology than 
Rl, R2 or R3 peaks which are due to the interaction between the guided or evanescent 
modes with the stub resonances (see previous Section). Oppositely, the LRl and LR2 
peaks correspond to crossing points in the band structure shown in Figure Eb. 

Very recent works [53j have shown several kinds of crossing points in anisotropic 
phononic crystals. As a general rule for these systems, the bands involved in the crossing 
point, or level repulsion state, exchange their polarization, so that the polarization 
remains a continuous function of the wave vector k. For the case of SC in which 
only longitudinal polarization is allowed, one can observe symmetric and antisymmetric 
acoustic fields. If the SC is symmetric with respect to the incident acoustic field, only 
symmetric acoustic fields will be supported. However, if the system is non symmetric, 
then both symmetric and antisymmetric fields can be supported. The stubbed SCW is 
an example of a non symmetric system. 

In our case, the points LRl and LR2 correspond to level repulsion states in the 
band structures that introduce a local band gap. The level repulsion state for these 
cases appears due to the crossing between an antisymmetric and a symmetric bands in 
a similar way as in phononic crystals thin plates [511 |52l |59] . Figures 16}:; and|6li represent 
detailed views of the band structures calculated using EPWE for the points LRl and 
LR2 respectively. One can see that, in both cases, the repelled bands are connected 
by means of an additional band, which is evanescent. Thus, EPWE does not predict a 
repulsion level state but an evanescent connection between bands. 

In order to analyze in more details the phenomenology involved in these level 
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Figure 9. (Color online) Analysis of the symmetry transfer between the repelled bands 
in stubbed waveguides. The analysis in (a) and (b) correspond to the repulsion states 
shown in Figure [6f and [6]2 respectively. The acoustic fields are represented using the 
absolute value of the pressure obtained from the Fourier transform of the eigenvectors 
for the eigenvalues marked with red squares. 



repulsion states, we have analyzed several modes for different frequencies all along the 
bands involved in the repulsion. Figure M shows the acoustic fields obtained from the 
Fourier transform of the eigenvectors for several modes in the repelled bands. The 
analyzed modes are marked with red squares which are connected using red arrows with 
their corresponding acoustic fields. 

One can observe that the modes in the repelled bands present a completely real 
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eigenvalue, k, whereas the connection bands have complex eigenvalue. The repelled 
bands are predicted by EPWE and they are exactly the same as the ones predicted 
by the PWE, however the difference appears in the evanescent connection between the 
repelled bands. 

First, we fix our attention on the acoustic fields in the repelled bands. One can 
observe in Figure [9] that for each band, depending on the position with respect to the 
evanescent connection, different symmetries of the field can be excited. For example 
points A and B that belong to the same band have different symmetries. A is on the left 
of the evanescent connection and B is on the right of the evanescent connection. The 
acoustic field corresponding to the point A is mainly sjTiimetric. However the acoustic 
field for the point B is completely antisymmetric. Then, inside the band, there is an 
exchange of the symmetry, so that the symmetry remains a continuous function of the 
wave vector k for each band. The same situations appear for the couple of modes D-E, 
F-G, an I-J. 

The repelled bands are connected by an evanescent band which is the responsible 
of the exchange of symmetry between bands. The acoustic fields of the modes inside the 
evanescent connection have a mixed symmetry. We have studied the modes C and H in 
the evanescent connections. First of all one can see the mixed symmetry of the acoustic 
field for the modes C and H. Moreover, one can observe the evanescent behaviour of 
the acoustic field inside the guide with a low decay rate indicating the small value of 
the imaginary part of the eigenvalue, k. It is worth to noting that the acoustic field 
predicted by MST for the frequencies corresponding to the modes C and H is fairly close 
to the one predicted by EPWE (Figure |8]). The MST predictions can be seen in the 
Figure [6^. 

4-2.1. Symmetric and antisymmetric stubbed SCW. Following the previous results, the 
repulsion states in SC should appear in systems presenting some degree of antisymmetry, 
as for example a stubbed SCW. In the benefit of the clarity of the previous results, we 
have analyzed two additional systems, a symmetric and another one with a degree of 
antisymmetry. Figure [TOb shows the complex band structure for a completely symmetric 
system with a symmetric stub as shown in the inset. Figure ITOb shows the complex 
band structure for the SCW with an antisymmetric stub which supercell is also shown 
in the inset. 

One can observe that the band structure of the symmetric stubbed SCW presents a 
crossing point but the bands are not connected, i.e., they are independent and only the 
symmetric one will contribute to the transmission properties of the system. However, 
for the case of the antisymmetric stubbed SCW one can observe again the evanescent 
connection between the repelled bands (blue diamonds). 
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Figure 10. (Color online) Complex band structures of two stubbed SCW, one 
symmetric (a) and the other one antisymmetric (b). The supercells for the EPWE 
calculations are shown in the insets. The evanescent connection between the repelled 
bands is marked with a blue diamonds for the case of the antisymmetric stubbed SCW. 



5. Concluding remarks 

The complex band structures obtained using tlie EPWE {k{u) metliod) sliow additional 
bands never revealed by the classical u{k) methods. In this work we have shown 
both theoretically, with two independent theoretical techniques (EPWE and MST), 
and experimentally, making use of SDReAMS, the interpretation of these evanescent 
additional bands. Due to the conservation of the overall number of bands for a 
determined frequency, interesting interpretations of the deaf bands and the level 
repulsion states have been obtained depending on the symmetry or not of the system. In 
the ranges of frequencies where a deaf band is traditionally predicted an evanescent mode 
with the excitable symmetry appears changing drastically the transmission properties 
of the system. On the other hand, we have interpreted, without loss of generality, 
the level repulsion between symmetric and antisymmetric bands in antisymmetric sonic 
crystals stubbed waveguides as the presence of an evanescent mode connecting both 
bands. These evanescent modes explain both the attenuation produced in this range 
of frequencies and the transfer of symmetry from one band to the other. Thus, the 
additional bands of EPWE are not an artifact of the calculation but a real physical 
behaviour can be attributed to them and as a consequence, the evanescent modes should 
be considered in the design of the systems based on periodicity. 
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